A NARX Model-Based Condition Monitoring Method for Rotor Systems

In this study, we developed a data-driven frequency domain analysis method for rotor systems using the NARX (Nonlinear Auto-Regressive with eXternal input) model established by system vibration signals. We propose a model-based index of fault features calculated in a multi-frequency range to facilitate condition monitoring of rotor systems. Four steps are included in the proposed method. Firstly, displacement vibration signals are collected at multiple monitored rotating speeds. Secondly, the collected signals are processed as output data and the corresponding input data is generated. Then, NARX models are developed with input and output data to characterize the rotor system. Finally, the NRSF (Nonlinear Response Spectrum Function)-based nonlinear fault index is calculated and compared to the healthy condition. An experimental application to the misaligned rotor system is also demonstrated to verify its effectiveness. Our results indicate that the value of the index directly reflects the severity of the misaligned fault.


Introduction
Rotating machinery, such as engines, fans, water pumps, etc., plays an important role in aerospace, energy, transportation, petrochemical, and other industrial fields [1][2][3]. With the increasing use of rotating machinery in mechanical equipment, researchers have become more concerned about its performance requirements. They strive to improve the machines' efficiency, reliability, and overall performance. However, the operating system of rotating machinery is often affected by various types of fault. Severe accidents can happen, resulting in economic losses. Therefore, it is necessary to monitor the status of the rotor system and detect early warning signs to avoid fatal damage.
Dynamic behavior can reflect the inherent condition of the structure system. Dynamic response changes can often be used as the point of structural damage detection. Previous studies showed that faults often introduce nonlinear dynamic characteristics to rotor systems; for example, the friction fault [4][5][6], misalignment fault [7][8][9], and crack fault [10,11] cause a multi-harmonic system output response other than the fundamental harmonic component. For this reason, the fault diagnosis method for rotor systems is often accompanied by a second harmonic vibration response. The more serious the misalignment fault is, the larger the proportion of the second harmonic [8,9]. Volterra series-based methods, such as GFRF (Generalized Frequency Response Function) and NOFRF (Nonlinear Output Frequency Response Function) [12], have also been used for the fault diagnosis of rotor systems [4,13,14]. This fault diagnosis method differs from observing the second harmonic, which only considers the output signal. NOFRFs or GFRFs are identified using the system outputs and inputs to ensure that the fault features are derived from the damage to the system rather than input excitation [4,15]. In addition, we propose a modified harmonic balance-alternating frequency/time domain method to analyze the nonlinear dynamics of a rotor system [16].
NOFRFs have been proposed as an alternative frequency domain representation of GFRFs, which circles GFRFs' multi-dimensionality problem. They are often difficult to measure and display in practice. Current approaches to Volterra series-based fault diagnosis of rotor systems also have some problems. NOFRFs require multiple tests with different amplitudes at the same speed [13]. Furthermore, data-driven-based NOFRFs evaluated from the NARX model, created using input and output data from the system, are only applicable to conditions where the input and output data are random signals [15]. It possesses the abundant frequency and amplitude characteristics of a NARX model, which can test the different properties of a system except for rotor systems. On the other hand, the truncation order of NOFRFs depends on their presentation ability and the extent to which NOFRFs can distinguish the structural properties of different damage conditions. In other words, although the NOFRFs up to the fourth order were considered in the present research, errors can still be caused if it is not appropriate a priori [14].
To solve the problems of a NOFRF fault diagnosis method, we developed a NARX model for rotor systems from the multi-harmonic signal generated by a speed-up process. We also used NRSF (Nonlinear Response Spectrum Function) [17], a system frequency analysis tool based on the notion that the system frequency response function consists of a first-order frequency response function and other orders combined. Then, we developed a condition monitoring method through an index calculated by NRSF to quantify the severity of the nonlinear fault. The proposed method is more efficient than observing the second harmonic method in the following experiment.
In summary, there are two innovations in this paper with respect to state of the art technology: (1) using Time Synchronous Averaging (TSA) technology to improve the quality of input and output signals; (2) proposing NFIs (Nonlinear Fault Indicators) to quantify the severity of the fault. This paper is organized as follows: Section 2 describes the theoretical background of the NARX model and the NRSFs method. The condition monitoring method for rotor systems is proposed in Section 3. Section 4 verifies the validity and feasibility of this method through an experiment. Finally, our conclusions are summarized in Section 5.

The NARX Model of a Nonlinear System
Compared to various data-driven models, NARX models are suited to conditionmonitoring rotor systems due to their physically interpretable nature based on being mapped into a frequency function [18]. The NARX model represents a large class of input-output nonlinear dynamic systems [19].
where F[•] is a nonlinear function that must be identified from the input u(k) and the output y(k); k represents the discrete time; n u and n y are the maximum input and output lags; and e(k) is the model prediction error assumed to be white noise.
Considering u(k) as a polynomial function, the NARX model (1) can be written as: where θ i 1 ···i l are the model coefficients, l being the degree of polynomial nonlinearity, n = n u + n y , and The polynomial NARX model (2) can also be represented as a linear-in-the-parameters model: where n M is the number of NARX model terms with full structure and n M = (n + L)!/(n!L!); p i (k), i = 1, . . . , n M are the monomials of the NARX model. Formula (2) is formulated by x 1 (k), . . . , x n (k); θ i , i = 1, . . . , n M . n M of the NARX model candidate is usually a high number. Identifying the NARX model helps select significant model terms from a large set of candidates based on input/output data so that a parsimony model structure can represent system dynamics. For example, a system with a nonlinear differential equation model: where u(t) and y(t) are the input and output; t is continuous time; and m, k 1 , c, and k 3 are the mass, linear stiffness, damping, and nonlinear stiffness, respectively. Formula (4) can be discretized into a difference equation through a forward-backward difference approach, where . y(t) and .. y(t) can be discretized as follows [20]: .
In practice, the NARX model identification typically uses the FROLS (Forward Regression Orthogonal Least Squares) algorithm by computing the ERR (Error Reduction Ratio) of each of the candidate terms [21], which will be introduced below.
The NARX model (4) can be written into a matrix form as:  where W = PA −1 , g = [g 1 , . . . , g n M ] T is an auxiliary parameter vector. Then the ERR criterion is defined as The significant model term with the highest ERR value is selected by comparing the ERR values of each column vector of the regression matrix P. The process is repeated to determine the n M significant model terms until where ρ is the threshold. In addition, an identified NARX model can easily transform into the frequency response function, revealing core invariant behavior and facilitating the analysis, design, and fault diagnosis for the system [15].

The NRSFs Representation of Nonlinear System
When a discrete-time NARX model is identified, the NRSF 1 (jω) [17], also known as the linear transfer function H 1 (jω) of the NARX model, can be directly evaluated through harmonic probing method [22] in the discrete-time domain.
where N u and N y represent lags in the linear model terms.
The Volterra functional series is shown in [18] y(t) = +∞ ∑ n=1 y n (t) (16) where y n (t) is the nth order output of the system. Fourier transform Y(jω) of y(t) can be expressed as: where ω is the frequency and Y n (jω) is the nth-order output frequency spectra obtained using Fourier transforms y n (t). Therefore, where, are defined as the NRSFs of nonlinear systems and NRSF ∞ (jω) represents the system of nonlinear dynamics [17]. Based on the NARX model (8), NRSF 1 (jω) can be evaluated as: Considering ω = [0, 20] rad/s, Figure 1 shows the NRSF 1 (jω) and NRSF ∞ (jω) of the NARX model under normalized amplitude excitation in the frequency range.   System damage usually introduces nonlinear dynamic characteristics in the structure, causing nonlinear behaviors in a system that originally presented linear characteristics [23,24]. For example, when rub-impact happens to rotor systems, the vibration waveform in the time domain is distorted because of the superposition of higher frequency components on the fundamental harmonic [4,5,7]. The vibration response of misalignment faults mainly consists of the fundamental harmonic and the second harmonic. The more severe the misalignment, the larger the proportion of the second harmonic [7,9,10]. The crack fault in the rotor system can also cause fundamental, second, and third harmonic components [10,11].
When an alternative system frequency domain expression differs from GFRF and NOFRF, NRSF divides the system's frequency response into linear and nonlinear parts without considering the GFRF dimension or the truncation order of NOFRF. The calculation results can be accurate and reliable with a higher calculation efficiency. Previous studies verified the sensitivity of NRSFs to changes in system parameters and system faults [17], whose results show that NRSFs can reflect the changes in system characteristics and provide effective feature extraction for system analysis, parameter estimation, and fault diagnosis.

Problem Statement
The condition monitoring method in this paper can be summarized into four steps: signal collecting, data processing, data-driven modeling, and fault feature extraction. The NARX model identification and NRSF method mentioned in Section 2 can be used to represent the rotor system and feature extraction of the model. Figure 2 shows the basic idea behind the proposed method. By extracting nonlinear characteristics from the NARX model established by the FROLS algorithm based on the vibration signals, which represent changes in physical characteristics caused by failure or damage in the rotor system, and comparing the fault index with the healthy condition, we can determine system faults System damage usually introduces nonlinear dynamic characteristics in the structure, causing nonlinear behaviors in a system that originally presented linear characteristics [23,24]. For example, when rub-impact happens to rotor systems, the vibration waveform in the time domain is distorted because of the superposition of higher frequency components on the fundamental harmonic [4,5,7]. The vibration response of misalignment faults mainly consists of the fundamental harmonic and the second harmonic. The more severe the misalignment, the larger the proportion of the second harmonic [7,9,10]. The crack fault in the rotor system can also cause fundamental, second, and third harmonic components [10,11].
When an alternative system frequency domain expression differs from GFRF and NOFRF, NRSF divides the system's frequency response into linear and nonlinear parts without considering the GFRF dimension or the truncation order of NOFRF. The calculation results can be accurate and reliable with a higher calculation efficiency. Previous studies verified the sensitivity of NRSFs to changes in system parameters and system faults [17], whose results show that NRSFs can reflect the changes in system characteristics and provide effective feature extraction for system analysis, parameter estimation, and fault diagnosis.

Problem Statement
The condition monitoring method in this paper can be summarized into four steps: signal collecting, data processing, data-driven modeling, and fault feature extraction. The NARX model identification and NRSF method mentioned in Section 2 can be used to represent the rotor system and feature extraction of the model. Figure 2 shows the basic idea behind the proposed method. By extracting nonlinear characteristics from the NARX model established by the FROLS algorithm based on the vibration signals, which represent changes in physical characteristics caused by failure or damage in the rotor system, and comparing the fault index with the healthy condition, we can determine system faults and their severity. However, there are still some related problems that must be solved before they can be used for condition monitoring of rotor systems in practical engineering.
Due to noise, the quality of the input and output signals significantly affects the NARX model's ability to characterize the system. Therefore, the measured signals should be processed before being used to establish a model.
The basic idea of most fault diagnosis methods is to extract the fault features and compare them with the healthy state to judge the current system condition [15]. Therefore, a fault indicator can help conduct the above process more conveniently.

Solution of Existing Problems
The random input that a NARX model often applies does not apply to rotor systems, which often operate at several rotating speeds [1][2][3]. The NARX model established using single harmonic signals can only be investigated under the same frequency [25]. Many researchers have analyzed rotor-bearing systems based on the speed-up process [21] and the frequency sweep system identification approach [26]. In this method, the NARX model of rotor systems is established using vibration signals under different rotating speeds, which reflect the frequency characteristics of the rotor system within the speed range.
If the rotating speeds of rotor systems in engineering practice are limited, then the range of monitoring speed Many data processing methods can address the noise in NARX modeling and enhance the NARX model's generalization ability, including regularization [27], AIC (Akaike Information Criterion), BIC (Bayesian Information Criterion), etc. [28]. Time Synchronous Averaging (TSA) technology [29] efficiently alleviates the influence of noise and measurement error on signals since the rotor system's vibration response signals are periodic. Before signal division, the vibration response signals 1 ( ) x t can be resampled with a sampling frequency higher than the original sampling frequency, making it more accurate to divide the signal into segments with more complete periods. The data processing method using TSA technology is presented as follows.  Due to noise, the quality of the input and output signals significantly affects the NARX model's ability to characterize the system. Therefore, the measured signals should be processed before being used to establish a model.
The basic idea of most fault diagnosis methods is to extract the fault features and compare them with the healthy state to judge the current system condition [15]. Therefore, a fault indicator can help conduct the above process more conveniently.

Solution of Existing Problems
The random input that a NARX model often applies does not apply to rotor systems, which often operate at several rotating speeds [1][2][3]. The NARX model established using single harmonic signals can only be investigated under the same frequency [25]. Many researchers have analyzed rotor-bearing systems based on the speed-up process [21] and the frequency sweep system identification approach [26]. In this method, the NARX model of rotor systems is established using vibration signals under different rotating speeds, which reflect the frequency characteristics of the rotor system within the speed range.
If the rotating speeds of rotor systems in engineering practice are limited, then the range of monitoring speed [ω min , ω max ] should be selected when the collecting is convenient. Collect the displacement vibration signal x 1 (t) under multiple monitoring speeds within [ω min , ω max ] at the same observation point, where i represents the ith monitoring rotating speed. Keep the observation point, collecting time t, and sampling frequency f s unchanged for each monitoring speed.
Many data processing methods can address the noise in NARX modeling and enhance the NARX model's generalization ability, including regularization [27], AIC (Akaike Information Criterion), BIC (Bayesian Information Criterion), etc. [28]. Time Synchronous Averaging (TSA) technology [29] efficiently alleviates the influence of noise and measurement error on signals since the rotor system's vibration response signals are periodic. Before signal division, the vibration response signals x 1 (t) can be resampled with a sampling frequency higher than the original sampling frequency, making it more accurate to divide the signal into segments with more complete periods. The data processing method using TSA technology is presented as follows.
Determine the number of periods N for each segment of data used for stacking and divide the signal into [M/N] segments, where M represents the total number of periods. Choose zero points at the start and end of each segment. Determine the sum of each segment of data and calculate the average as follows: where y i (t) is the processed data of the ith monitoring speed, x i,n (t) represents the nth segment of data for the ith monitoring speed. Considering the duffing system (5), the output response under u(t) = 2 sin(ω 1 · t), where ω 1 = 15rad/s, demonstrates the signal process method. To highlight improvements in the processed signal, a noise e(t) is mixed with the system output x 1 (t) by setting e(t) as a uniformly distributed random sequence in [−4, 4] × 10 −4 . We collected 300s of the system output with a sampling frequency of f s = 1024 and obtained M = 714 periods to conduct the TSA method. Figure 3 shows the number of periods for each segment is N = 5 where an interval of two stars denotes a period. The signal is divided into [M/N] = 142 segments. After summation averaging (15), Figure 4 shows the output y 1 (t) of using TSA technology. in the processed signal, a noise ( ) e t is mixed with the system output 1 ( ) x t by setting ( ) e t as a uniformly distributed random sequence in 4 [ 4,4] Figure 3 shows the number of periods for each segment is 5 N  where an interval of two stars denotes a period. The signal is divided into segments. After summation averaging (15), Figure 4 shows the output 1 ( ) y t of using TSA technology.
One segment  The rotor system is affected by a centrifugal force from an unbalanced mass. However, the centrifugal mass and distance are sometimes unknown or difficult to measure. Therefore, the input ( ) i u t for the NARX model can be generated in a harmonic form with a normalized amplitude.
NARX models can adapt to different system parameters by changing the coefficients of model terms so that the amplitude of the input will not affect its ability to represent the rotor system. At each rotating speed, the sampling frequency, vibration frequency, and data length of the input signal are consistent with the output signal, and the starting point is also set to zero. The sampling frequency can be appropriately reduced after data processing to create NARX models more efficiently. ( ) e t as a uniformly distributed random sequence in [ 4,4] Figure 3 shows the number of periods for each segme 5 N  where an interval of two stars denotes a period. The signal is divided [ / ] 142 M N  segments. After summation averaging (15), Figure 4 shows the ou 1 ( ) y t of using TSA technology.  The rotor system is affected by a centrifugal force from an unbalanced mass. H ever, the centrifugal mass and distance are sometimes unknown or difficult to mea Therefore, the input ( ) i u t for the NARX model can be generated in a harmonic form a normalized amplitude.

One segment
NARX models can adapt to different system parameters by changing the coeffic of model terms so that the amplitude of the input will not affect its ability to represen rotor system. At each rotating speed, the sampling frequency, vibration frequency, data length of the input signal are consistent with the output signal, and the starting p is also set to zero. The sampling frequency can be appropriately reduced after data cessing to create NARX models more efficiently. The rotor system is affected by a centrifugal force from an unbalanced mass. However, the centrifugal mass and distance are sometimes unknown or difficult to measure. Therefore, the input u i (t) for the NARX model can be generated in a harmonic form with a normalized amplitude.
NARX models can adapt to different system parameters by changing the coefficients of model terms so that the amplitude of the input will not affect its ability to represent the rotor system. At each rotating speed, the sampling frequency, vibration frequency, and data length of the input signal are consistent with the output signal, and the starting point is also set to zero. The sampling frequency can be appropriately reduced after data processing to create NARX models more efficiently.
In this case, several data sets under different monitoring speeds must be handled simultaneously and one common NARX model is required. This can be achieved by the data splicing method, where the inputs u i (t) and outputs y i (t) under different speeds are spliced into U(t) and Y(t), respectively, in the same order. Then, the FROLS algorithm can be applied to establish NARX models. To quantify the intensity of nonlinear faults and facilitate judgment of the rotor system condition, NFIs (Nonlinear Fault Indicator) are proposed in this method, which are defined as follows: (23) where ω min and ω max are the minimum and maximum values of the frequency range, respectively. It is worth noting that NFI is a special case of NOFRF-based indices [30].
Under the investigation frequency ω n selected in [ω min , ω max ], we applied a single harmonic excitation with the normalized amplitude sin(ω n t) to the NARX model to calculate NRSF.
Compared to the fault features calculated at a single frequency, the index proposed in this method contains more information [15], resulting in a more reliable monitoring result.

Condition Monitoring Procedure
The condition monitoring method for rotor systems using NARX model and NRSFbased analysis can be summarized as follows: Step 1: Collect the displacement vibration signals x i (t) under multiple monitoring rotating speeds ω i , where i represents the ith monitoring rotating speed; Step 2: Process the collected signals using the TSA method and generate the corresponding input signal. Then, splice u i (t) and y i (t) into U(t) and Y(t), respectively; Step 3: Establish a NARX model to characterize the rotor system based on input U(t) and output Y(t) data using the FROLS algorithm; Step 4: Calculate the nonlinear fault index (17) based on the identified NARX model and compare the results to healthy conditions.

The Misaligned Rotor System Experiment
We conducted this experiment to demonstrate the condition monitoring method using an experimental rotor system with various levels of misalignment ( Figure 5). Step 2: Process the collected signals using the TSA and data splicing methods.
LMS mobile front end  By setting different offset distances between the edge of the bearing seat and the pedestal, we achieved different levels of misalignment fault. Each increase of 1 mm offset distance was set as a misalignment level, and a total of four levels were set in the experiment, namely level 1, level 2, level 3, and level 4, respectively. The eddy current displacement sensor was aligned with the horizontal direction of the rotating shaft and connected to the LMS mobile front end to collect the experimental rotor system's response.
Using the proposed method in Section 3.3, we performed condition monitoring of the rotor system under the level 1 misalignment fault.
Step 1: Collect the displacement vibration signals. The displacement vibration signals were collected under rotor system speeds of ω 1 = 1550 r/min, ω 2 = 1650 r/min, and ω 3 = 1700 r/min at the same observation point, denoted as x 1 (t), x 2 (t), and x 3 (t), respectively, with a sampling frequency f s1 = 2048 Hz and collecting time of t = 33 s; Step 2: Process the collected signals using the TSA and data splicing methods. The vibration response signal under ω 3 = 1700 r/min at the level 1 misalignment demonstrates data processing in this method. The collected vibration response signal x 3 (t) was interpolated by a cubic spline and sampled again with the new sampling frequency f s2 = 20480 Hz, which was L = 10 times the original sampling frequency f s1 . The number of periods for each segment was determined to be N = 20 and the whole M = 895 periods were divided into [M/N] = 44 segments. Each segment of data was summed together and an average was calculated. After the system outputs were processed, the inputs for NARX modeling were generated using (16).
The data splicing method was applied to handle the three data sets. Under the level 1 misalignment of the rotor system, the input and output data were sequentially spliced, and the data after splicing was resampled with the sampling frequency f s3 = 1024 Hz. Figures 6 and 7 show the splicing results of input U(t) and output Y(t), respectively.  Step 3: NARX modeling using the FROLS algorithm. By using the FROLS algorithm, we initialized the maximum time delays for the system input n u = 3, output n y = 3, and the degree of polynomial nonlinearity l = 3. Based on the input U(t) and output Y(t) data after splicing, the NARX model for level 1 of the misaligned rotor system is: Step 3: NARX modeling using the FROLS algorithm. By using the FROLS algorithm, we initialized the maximum time delays for the system input 3 u n  , output 3 y n  , and the degree of polynomial nonlinearity      Figure 9 shows NFI (17) calculations for the NARX model of the rotor system in the first misalignment level. Using the same method as steps 1-4, we obtained the NFI of the rotor system in levels 2-4 of misalignment fault ( Figure 8). As shown in Figure 9, the NFI values at four misalignment levels are 0.076, 0.21, 0.25, and 0.45, respectively. These values have an increasing monotonic trend with the severity of the misaligned fault, meaning that the index NFI's value can directly reflect the severity of the misaligned fault. Therefore, the method proposed in this paper accurately identifies rotor system faults. Step 4: Calculate the NRSF-based nonlinear fault index. The frequency domain characteristics of the NARX model were investigated in the range of monitoring rotating speeds [1550, 1700] r/min. The investigation frequencies ω n used to calculate NRSF were set at every other 0.5 rad/s within the range [1550, 1700] r/min. Under the investigation frequency, a single harmonic excitation was applied to NARX models to calculate NRSF. Figure 9 shows NFI (17) calculations for the NARX model of the rotor system in the first misalignment level.   Figure 9 shows NFI (17) calculations for the NARX model of the rotor system in the first misalignment level. Using the same method as steps 1-4, we obtained the NFI of the rotor system in levels 2-4 of misalignment fault ( Figure 8). As shown in Figure 9, the NFI values at four misalignment levels are 0.076, 0.21, 0.25, and 0.45, respectively. These values have an increasing monotonic trend with the severity of the misaligned fault, meaning that the index NFI's value can directly reflect the severity of the misaligned fault. Therefore, the method proposed in this paper accurately identifies rotor system faults.

Comparison to the Traditional Method
Misalignment faults are often accompanied by second harmonic vibration responses. The more severe the misalignment fault is, the larger the proportion of the second har- Using the same method as steps 1-4, we obtained the NFI of the rotor system in levels 2-4 of misalignment fault ( Figure 8). As shown in Figure 9, the NFI values at four misalignment levels are 0.076, 0.21, 0.25, and 0.45, respectively. These values have an increasing monotonic trend with the severity of the misaligned fault, meaning that the index NFI's value can directly reflect the severity of the misaligned fault. Therefore, the method proposed in this paper accurately identifies rotor system faults.

Comparison to the Traditional Method
Misalignment faults are often accompanied by second harmonic vibration responses. The more severe the misalignment fault is, the larger the proportion of the second harmonic [8,9]. Therefore, the amplitude of the second harmonic was always used to characterize misalignment faults in rotor systems. We compared methods to verify the advantages of the method proposed in this paper.
The vibration response spectrums of each misaligned rotor system level were calculated using the experimental data investigated beforehand. Figure 10 shows the amplitude of the second harmonic. In the traditional method, there is no obvious upward monotonic trend with the severity of the misaligned fault. For example, the amplitude of the second harmonic is gradually increasing with the severity of the misaligned fault in levels 1-3 but decreases at level 4.  The diagnosis method based on output signals is affected by system input and noise and cannot reflect the essential characteristics of the system. By contrast, we propose a model-based fault feature extraction method with an index calculated in a multi-frequency range, which is more reliable for condition-monitoring rotor systems.

Conclusions
In this study, we propose a NARX model-based NRSR method to facilitate condition monitoring for rotor systems. We used TSA technology to improve the quality of training data for NARX modeling. We also used an NFI indicator to quantify fault severity in rotor systems. Our specified conclusions are as follows: (1) NRSF-based methods are more suitable for situations where the prior conditions are insufficient and the fault characteristics are unknown since they do not have to consider the truncation order of NOFRF; (2) The NFI values at four misalignment levels were 0.076, 0.21, 0.25, and 0.45, which have a monotonic trend with the severity of the misaligned fault validating the proposed method's efficiency over the second harmonic method.  The diagnosis method based on output signals is affected by system input and noise and cannot reflect the essential characteristics of the system. By contrast, we propose a model-based fault feature extraction method with an index calculated in a multi-frequency range, which is more reliable for condition-monitoring rotor systems.

Conclusions
In this study, we propose a NARX model-based NRSR method to facilitate condition monitoring for rotor systems. We used TSA technology to improve the quality of training data for NARX modeling. We also used an NFI indicator to quantify fault severity in rotor systems. Our specified conclusions are as follows: (1) NRSF-based methods are more suitable for situations where the prior conditions are insufficient and the fault characteristics are unknown since they do not have to consider the truncation order of NOFRF; (2) The NFI values at four misalignment levels were 0.076, 0.21, 0.25, and 0.45, which have a monotonic trend with the severity of the misaligned fault validating the proposed method's efficiency over the second harmonic method.